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Abstract 

We propose to measure the Cabibbo-Kobayashi-Maskawa parameter 2/3 + 7 using B° decays 
involving several intermediate states, and describe a general formalism that applies to a broad 
class of decays. The main advantage of this method is that the ratios between the interfering 
amplitudes can be measured without requiring external input. In addition, discrete ambiguities 
are resolved. 
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I. INTRODUCTION 



CP violation is one of the most important topics in current particle physics research. 
In the standard model, CP violation arises due to a single complex phase in the Cabibbo- 
Kobayashi-Maskawa matrix V [1]. A major goal of B meson physics is to measure the angles 
and sides of the CKM unitarity triangle. Theoretically clean measurement methods are 
crucial for obtaining these parameters accurately. The BaBar [2] and Belle [3] measurements 
of the parameter sin(2/3), where (3 = arg (— V c dV* b /V t dV t * b ), confirm the standard model to 
within the precision of the experiments, and increased precision is expected in the coming 
years. 

Crucial studies of the CKM mechanism and constraints on new physics can be obtained 
by measuring the CKM angle 7 = arg (— V u dV* b /V c dV* b ). The greatest challenges presented 
by these measurements is that they require very large data samples and are subject to 
discrete ambiguities. It is therefore important to use every possible mode and method for 
measuring 7, and to devise methods that help resolve the ambiguities. 

An important class of measurements makes use of decays such as B — > D~n + to measure 
2/3 + 7. Proposed initially by Dunietz [4], the first attempts to measure time-dependent CP 
asymmetries proportional to sin(2/3 + 7) and cos(2/3 + 7) have been conducted by BaBar [5] 
and Belle [6] using the modes B — > D^~it + and B — > D~p + . While these measurements 
are currently statistically limited, their precision will become significant as more data are 
accumulated. At that stage, the greatest difficulty in extracting 2/5 + 7 from these results 
will be the lack of precise knowledge of the ratio between the interfering amplitudes, defined 
as r = \A(B° -> D^h + )/A(B° -> D^~h + )\, where h+ indicates the light hadron tt+ or 
P + - 

In principle, r may be obtained from the difference between the magnitudes of two terms 
with different time dependences in the decay rate. The relevant terms are (1 + r 2 ) and (1 — 
r 2 ) cos Amt, where Am is the B° — 5° oscillation frequency. However, with r ~ 0(1 — 2%), 
extracting it from the 0(r 2 ) difference between these 0(1) terms requires prohibitively large 
data sets. Thus, the time-dependent measurement has negligible sensitivity to the value of 
r, which must therefore be obtained by assuming factorization and SU(3) symmetry to 
make use of the ratio of branching fractions B(B° — > h + )/B(B° — * D ( -*^~h + ). This 
approximation ignores the contribution of annihilation diagrams and some SU(3) breaking 
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effects, and is taken to have a theoretical error of roughly 30% [5]. 

In B — > D*~p + , the single parameter r is replaced by a matrix p mn of ratios between 
the magnitudes of the b — > ucd and 6 — > cud contributions of the of three different helicity 
amplitudes contributing to the decay. It has been shown [7] that p mn may be obtained 
using only first-order 0(p mn ) terms. Not having to rely on small second-order terms or 
external input regarding amplitude ratios, this provides a much improved, theoretically 
clean measurement of 2/3 + 7. 

In this paper we generalize and extend that method to other decays that proceed through 
more than one intermediate state. Examples include B — > D~p + , which can interfere with 
B — > D~p + (1450) and non-resonant B — » D~7t + tt°; B — > D^a^ , where non-resonant con- 
tributions are expected under the a\ peak; and B — > D**~tt + , where interference between 
several excited charmed mesons may be realized in the decays D**~ — > Dir and D**~ — > D*n, 
in addition to possible contributions from non-resonant decays. 

In all these cases, the interfering contributions have overlapping yet different distributions 
in relevant analysis variables. The first of these variables is the invariant mass squared s 
of the final state of the resonance. The second variable s' typically describes an angular 
distribution that is fully determined by the spin of the resonance. In the case of B — > D~a^ , 
s' corresponds to the two variables of the Dalitz plot of the a\ decay. 

Our method applies equally well to modes with higher excitations, such as B — ► D*~p + , 
B — > D*~af , B — » D**~ p + , and B — > D**~a^ , in which s' corresponds to several angular and 
mass-related variables. In addition to the interference between several resonances and non- 
resonant contributions, these decays involve several helicity amplitudes, which are treated 
as different intermediate states in our method. 

II. MEASURING 2/3 + 7 

Let us consider a decay of the type described above, involving the interference of N 
intermediate states. We denote the final state by / (/) if it contains a c (c) quark. The four 
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decay amplitudes of interest are 

N 

A(B° - /) = A(B° - 7) = £ ^')e iAra , 

m=l 

iV 

m=l 

N 

A(B° - /) = £ a m ^ m ( S , S ')e i(5 ™- 7) , (1) 

m=l 

where A m (5 m ) is the CP-conserving phase and v4 m (a m ) is the magnitudes of the b — > 
cwci (6 — > ucd) decay amplitude proceeding via intermediate state m, and g m (s, s') is a 
known function of the final state variables s and s' that depends on the nature of the 
intermediate state m. For example, for / = D~7t + tt° and m being the index of the D~ p + 
intermediate state, g m (s, s') = R(s)s', where s is the square of the im° invariant mass, R(s) 
is a Breit-Wigner function, and s' is the cosine of the angle between the momenta of the B 
and of one of the pions, calculated in the ir + ir rest frame (the "helicity" angle). Vector- 
vector intermediate states, such as D*~p + , must be further divided into the different helicity 
amplitude, each of which has a different s' dependence. 

With the above equations, the time- dependent decay rates for B°(t) — > f and B°(t) — > f 
become 

T(B°(t) - /) = e~ rt [tmn + C mn cos(Amt) - S~ n sin(Amt)] , (2a) 

m,n 

T(B°(t) - /) = e~ rt Vmn - C mn cos(Amt) - S+ n sin(Amt)] , (2b) 

m,n 

where for convenience we define the symbols 

1 '(A A ) '(S <5 ) 

Iran = ^ \_9m9n i^-^mAnC ~\~ Cl m Cl n G ) } ' 

C = -in n* ( A A P -*(A„-A m ) _ p -i(K~S m )\\_ 

°mn — 1111 \ymy n n -nU'm ( i & J i 

Q+ = Tm in n* A n p -*( 5 n-A TO ) i<H 

= -(2/3 + 7). 

(3) 

The decay rates for 5° decays are obtained from the B° rates by inverting the sign of 

4 



the cos(Amt) and sin(Amt) terms. They double the statistics but do not yield additional 
information. 

Next, we determine the conditions under which all the unknown parameters of Eqs. (2) 
can be obtained from the measurement, and show that these conditions are satisfied in the 
typical case of interfering Breit-Wigner resonances and a possible non-resonant contribution. 

The three terms of Eq. (2a) are distinguishable based on their different time dependences, 
thus determining their coefficients. The relative differences between I mn and C mn are of order 
(a m a n ) I (A m A n ) ~ r 2 ~ 10~ 4 , which is practically unobservable. As a result, these terms 
yield the parameters A m and A m , while a m and 6 m are measured from the coefficients of the 
sin(Amt) terms, as described later. To study the conditions for obtaining A m and A m , we 
expand 

EA A in n *p-^ An - Am H - V In I 2 A 2 
m.n m 

+ 2 ^ ^{g m g* n )A m A n cos(A n - A m ) 

+ 2 Im(gmg* n )A m A n sin(A n - A m ). (4) 

If \g m \ 2 , Re(g m g^), and lm(g m g^) all have different s and/or s' dependences, Eq. (4) yields 
N 2 unique observables, which is more than enough to determine the 2N — 1 unknowns A m 
and A m (one of the A m phases is a global phase and can be chosen arbitrarily) for N > 2. 
This uniqueness condition is satisfied when all the g m are Breit-Wigner functions, 

9m[S) = s-M 2 m + l M m Tj (5) 
even when all contributions have the same s' dependence. A non-resonant contribution 
gi = 1 introduces N — 1 degenerate relations: 

lm (9i9 m ) = \9m\ 2 , (6) 

where g m (m > 1) is a Breit-Wigner function. In this case, the number of observables is 
reduced to iV 2 — (N — 1). However, a solution still exists for N > 2, and this solution is 
unambiguous when the non-resonant contribution is small enough relative to the resonant 
contributions. In addition, most practical cases involve resonances with total spins differ- 
ent from 0, and hence s' dependences that distinguish them from a non-resonant s-wave 
contribution. This guarantees a unique solution of Eq. (4) in terms of A m and A m . 
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We note that these conclusions do not depend on the assumption that the a m a n terms 
in I mn and C mn are negligible. In fact, they apply equally well to the a m a n {g m g^ l e~ % ^ n ~ &mS> } 
terms in Eq. (3). 

We now show how the coefficients of the sin (Ami) terms in Eq. (2) yield the values of 
the remaining 2N + 1 unknowns, namely, a m , 8 m , and 0. The coefficients are 

^SZn = A m a m \g m \ 2 sin(0 ± S mm ) 

mn m 

in 

a n cos(0 ± 5 nm ) ± A n a m cos(0 ± 5 mn )] 

m<n 

+ ^ R e (9m9n) [A m a n sin(0 ± S nm ) + A n a m sin(0 ± 8 mn )\ (7) 

m<n 

where 5 nm = S n - A m . If \g m \ 2 , Re(g m g*), and lm(g m g*) are all different, Eq. (7) yields N 2 
observables for S~ n and N 2 for S+ n - It is therefore possible to obtain all the unknowns for 
N > 2. 

III. DISCRETE AMBIGUITIES 

In the N = 1 case, only the first line in Eq. (7) is non-vanishing. The measurement of <f) 
then suffers from an eight-fold ambiguity, due to the invariance of the observable sin(</>±5 mn ) 
under the three symmetry operations [8]: 

Sir/2 = 4> —> ^mm. + 7r/2 , 5 mm — > (j) - 7r/2, 

S'k = S mm + 71 , 5 mm -> + 7T, (8) 
S'± =0 > 71 — <f) , 5 mm > — S mm . 

In the typical N > 1 case, S' 7r /2 and S± are no longer good symmetries, since they are 
broken by the cos(0 ± 5 mn ) terms. Furthermore, these terms are distinguishable from the 
sin(0± 5 mn ) terms by virtue of the different s and/or s' dependences of Re(g m g*) and \g m \ 2 
or lm(g m g^). This further improves the measurement of 0. 

IV. DISCUSSION AND SUMMARY 

In this paper, we have outlined the formalism for measuring 2/9 + 7 with neutral B meson 
decays involving interference between several intermediate states. We have shown that, 
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despite involving a more complicated analysis, these decays have distinct advantages over 
B — > D^~n + , once our formalism is applied to their analysis, thus enhancing the overall 
precision with which 2(3 + 7 is known. 

First, as already noted for the special case of B — > D*~V + decays [7], our method is 
sensitive to 2(3 + 7 using only first order terms in the ratios a m /A n between the b — > ucd 
and b — > cud amplitudes. By contrast, in B — > D^~tt + , or in the analysis of other decay 
modes that ignores the contribution of multiple intermediate states, one needs to extract 
r = a\/A\ from 0(1 — r 2 ) terms, or rely on external measurements and incur a large 
theoretical uncertainty. Since r is as small as 1 — 2%, this advantage is realized in our method 
even when the amplitude of one of the interfering intermediate states is much greater than 
the others. 

Second, the B — > D*~n + measurement is subject to an eight-fold ambiguity, while in our 
method, the ambiguity is only two-fold. 

Our method is not completely model-independent, since one has to assume specific forms 
for the g m functions, such as a Breit-Wigner for the resonances. However, this model depen- 
dence is much smaller than the 30% theoretical error estimated for r. Most resonances are 
well understood, and their shapes can be studied with the terms of Eqs. (4). In addition, 
the number of observables in Eq.(7) is greater than the number of unknowns when there 
are more than two intermediate states. The additional constraints may be used to further 
reduce the model dependence associated with some g m functions. 

We emphasize that these conclusions and the formalism presented here do not depend on 
a specific final state, but apply whenever enough is known about the g m (s, s') functions for 
a solution to be obtainable, which in practice holds for a majority of the cases. 

It is interesting to note some similarities and differences between the method we present 
here and methods developed for measuring 7 in B — »• DK. Multi-body final state B — > DK 
decays (such as B~ -> DK~n° [9], B~ -> D**K~ [10], B~ -> D*K*~ [11], or B -> DK 
with multi-body D decays [12]) have been shown to improve the measurement of 7. This 
improvement is mostly due to the resolution of ambiguities and the ability to make efficient 
use of many B and D modes. As we have shown here, similar advantages are realized by 
interference between intermediate states in the measurement of 2(3 + 7 with multi-body 
B — > .D~7r + -like modes. But in addition, these measurements benefit mostly from the fact 
that they do not depend on the very small r 2 terms. By contrast, B — > DK decays are 
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governed by the amplitude ratio r B = \A(B~ -> D K~)/A{B- -»■ ~ 10 - 20%, 

which is about an order of magnitude larger than r. Therefore, the sensitivity advantage 
brought about by interference between intermediate states is much greater in B — > D~tt + 
than in B -> Dif . 
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